Evolution of Anisotropies in Eddington-Born-Infeld Cosmology 
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Recently a Born-Infeld action for dark energy and dark matter that uses additional afiine con- 
nections was proposed. At background level, it was shown that the new proposal can mimic the 
standard cosmological evolution. In Bianchi cosmologies, contrary to the scalar field approach (e.g., 
Chaplygin gas), the new approach leads to anisotropic pressure, raising the issues of stability of 
the isotropic solution under anisotropic perturbations and, being it stable, how the anisotropies 
evolve. In this work, the Eddington-Born-Infeld proposal is extended to a Bianchi type I scenario 
and residual post-inflationary anisotropies are shown to decay in time. Moreover, it is shown that 
the shears decay following a damped oscillatory pattern, instead of the standard exponential-like 
decay. Allowing for some fine tuning on the initial conditions, standard theoretical bounds on the 
shears can be avoided. 
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I. INTRODUCTION 

One of the greatest puzzles in current fundamental 
physics is the nature of dark energy and dark matter 
d 0]. Observationally their existence is well rooted in 
the results of many independent experiments, e.g. [1, H[, 
while theoretically the understanding of these entities is 
far from satisfactory. The simplest way to model the 
dark sector, i.e. dark energy and dark matter, in good 
agreement with observations is through the now stan- 
dard ACDM model 0. Following this approach, at back- 
ground level, dark matter enters as an additional dust 
contribution and dark energy as a positive cosmological 
constant. 

Considering their apparent dissimilarities, that dark 
matter and dark energy might be two faces of a single 
mysterious entity is an appealing proposal. This direc- 
tion was first explored through the employment of the 
Chaplygin gas and its generalized version[5|, which use a 
(generalized) Born-Infeld action for scalar fields. Albeit 
compatible with with cosmological data at background 
level [6], once density perturbations are considered the 
Chaplygin gas approach has no satisfactory answer to 
the structure formation and the cosmic microwave back- 
ground power spectrum Q . 

Recently, motivated from a general relativity extension 
to spaces whose space-time metric can be degenerated, 
a new model that describes an unified dark sector was 
proposed |8|. In this Eddington-Born-Infeld model (or 
Einstein-Eddington-Born-Infeld), the space-time metric 
interacts with the Eddington action [9( through a Born- 
Infeld coupling. 

The Eddington action was proposed as an alternative 
action for gravitation [9(. Its fundamental fields are the 
affine connections and no metric appears in the action. 



It reads 



In above, T 



5 Edd [r] = / J\R^(r)\d 4 x 



(i) 
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aj3 — r^ a are independent afiine connections 
and \R^ V \ is the modulus of the determinant of the sym- 
metric Riemann curvature tensor. The action (fT]) docs 
not depend on the space-time metric but it is equivalent 
to the Einstein-Hilbert action in the presence of a non- 
null cosmological constant [l(| ■ Considering that current 
cosmological observations favors a non-null cosmological 
constant, the latter point is welcome. Nevertheless, the 
coupling of gravity in this formulation to the rest of the 
universe is not a straightforward issue. The fundamental 
fields are the affine connections, but direct couplings of 
these to other structures are typically inconsistent, since 
the former are not tensors [271 ]. 

The above property turns to be a welcome feature in 
the Eddington-Born-Infeld model. In the latter, the ad- 
ditional affine connections that come from the Eddington 
part are interpreted as the dark sector fundamental fields, 
which couples to the rest of the universe only through 
gravity. 

Besides the aforementioned motivation, the employ- 
ment of additional connections in cosmology, instead of 
scalar fields, is a new curious possibility which deserves 
to be further investigated. The Eddington-Born-Infeld 
model can be seen as an afiine connection version for the 
Chaplygin gas approach, where the scalar kinetic term 
inside the determinant is replaced by the curvature ten- 
sor for the additional "dark" connections. 

The purpose of this work is to analyze the Eddington- 
Born-Infeld cosmology beyond the isotropic and homo- 
geneous background. Of particular interest is post- 
inflationary shear perturbations, which will be studied in 
a Bianchi I framework. It generalizes the flat Friedmann- 
Lemaitre-Roberton- Walker cosmology to homogeneous 
and anisotropic spaces without rotation. For reviews on 
Bianchi cosmologies see [HI. For some recent develop- 
ments on the Bianchi I phenomenology see [13 . [Tij and 



2 



references therein. While homogeneous scalar fields in 
Bianchi backgrounds exert isotropic pressure, preserving 
the shears tendency to decay, the approach with addi- 
tional affine connections generate anisotropic pressure for 
backgrounds not exactly isotropic. Hence, a priori, the 
latter approach could lead to an anisotropic universe in 
conflict with observations [H, EB|, even if the universe 
anisotropy just after inflation was negligible. In the fol- 
lowing sections it is shown that this a priori drawback for 
using additional connections in cosmology is not present 
in the Eddington-Born-Infeld model, the shears also de- 
cay in this model (although not monotonically) . In par- 
ticular this implies that the isotropic solution found in 
Q is stable under shear perturbations. Moreover, with 
the purpose of disclosing model predictions, the way the 
shears decay is also analyzed, and it is shown that it 
follows a characteristic pattern. 

In the following section we review the Eddington-Born- 
Infeld model, establish notation and introduce some pro- 
cedures that will be useful for the following sections. 
Sec. Ill is devoted to the study general properties of 
the model, for both late times and early times, that do 
not rely on the initial conditions that parametrize the 
anisotropics evolution. Sec IV presents numerical solu- 
tions for the Eddington-Born-Infeld model in a Bianchi I 
background, it illustrates Sec. Ill results, shows the evo- 
lution of anisotropies at small redshifts — a region that 
the Sec. Ill approach is inconclusive — and presents 
two peculiar solutions that can avoid standard theoreti- 
cal bounds on the shears [Hi [HI. Finally, in Sec. V, our 
conclusions and perspectives are presented. 



II. REVISITING THE 
EDDINGTON-BORN-INFELD ACTION AND ITS 
HOMOGENEOUS AND ISOTROPIC 
COSMOLOGY 

This section reviews some key results of [8[ and in- 
troduces procedures that will be useful to the following 
sections. 

The action proposed in reads 



S[g, C, *] = J £ m (*,g)d i x + 



(2) 



\g^\R- 



d x. 



where {g^v) is the space-time metric, k = l/(8nG), W 
stands for any additional matter fields, l-X^I is the ab- 
solute value of the (X^) determinant, is the sym- 
metric Ricci tensor constructed with the symmetric con- 
nection CP and R is the standard curvature scalar con- 
structed with the space-time metric and the metric con- 
nection T(g). The dimensions of the constants I 2 and a 
read [I 2 ] = m~ 2 , [a] = 1. The affine connection cou- 
ples to g^u through a Born-Infeld interaction, and only 



influence additional fields indirectly through the coupling 
of g^ u to them. 

In the absence of C m , the equations of motion of the 
action © read 



G ^ ~ al 2 



1 \g- l 2 K 



i 9tia {g-l 2 K)^gp v , (3) 



^^\g-l 2 K\(g~l 2 I<r") =0, 



(4) 



where the indices inside the determinants were omitted, 
((g — l 2 K) a/3 ) is the inverse of the matrix {g a p — l 2 K a p) 
and V c is the covariant derivative constructed with the 
C connection. 

The equations of motion of the pure Eddington ac- 
tion |T]) are given by @ without g^v- In order to ex- 
press such equations in a form closer to the Einstein 
equations, one introduces an auxiliary non-degerate sym- 
metric tensor q^ v such that V^q^ = 0, which implies 
that Cp p can be written as functions of q^ and deriva- 
tives. Letting q^ oc K^, Eq. (without g^) becomes 
an identity and the equations of motion are now simply 
Qfiv A = K/j, v (q), where A is a constant. Interpreting q M „ 
as the space-time metric, one sees that the latter equa- 
tions are just Einstein equations for a de Sitter space-time 
with A as the cosmological constant. For other details on 
the correspondence of this formulation with the standard 
one, see [10| . 

A similar trick is also useful for the action ([2]). The 
main differences are that the auxiliary tensor g M „ is not 
the space-time metric, and that the value of the arbitrary 
constant A corresponds to a rescaling of q^ v that has no 
dynamical role. The resulting new expressions for the 
equations of motion are 



Gnu = -A 



Q\ a/3 



(5) 



K^ v {q) =A[aq fJiU + (1 - a) g^]. (6) 

In above we have re-parametrized the model from {I 2 , a} 
to {A, a} using the same relation of [|[, 

/2 = T^A' W 

1 — a A 

but a different normalization for q^ v was selected. To 
change to the original normalization of @| one should 
replace every g M „ by (1 - a) q^. 

The simplest nontrivial solution of Eqs. (JHJ © are 
found for C? = This implies that q^ = jg^v, 

for constant 7. From (©, one concludes that this corre- 
sponds to a de Sitter solution with A7 as the cosmological 
constant; and requiring compatibility with ^ leads to 
7=1. The latter defines the selected q^ v normalization. 

Eq. ([6]) can be set in the standard form of Einstein 
equations, 



A(l-a) 



1 a/3 \ ( Q \ 

a — 1 I / 



Now we turn to the isotropic and homogeneous cos- 
mology spawned by the equations 0) with 

(g^)=dmg(-N 2 a 2 a 2 a 2 ), (9) 

(^)-diag(-X 2 Y 2 Y 2 Y 2 ). (10) 

All the above variables depend on time alone. The equa- 
tions of motion, in the presence of usual matter and ra- 
diation, can be stated as 



3H 2 = N 2 A 



■mi = 



NY 3 



Xa 3 
A(l - a)X 2 



Pb 
k 



k 



(11) 



2a N 2 a 2 . 

T^-^ + 3 ^)' (12) 



3a z H 



XY d /NY 3 \ 



N 



dt\X 



(13) 



where H = a/ a, Hy = Y/Y, pb is the baryon energy 
density and p r is the radiation energy density. Once the 
temporal gauge is fixed, one has three equations for three 
unknowns. The above equations are the Goo (fTT|) con- 
straint, the Qoo constraint (fT2|) . and the energy momen- 
tum conservation for the "dark fluid" (TIB")) . The latter 
can be found from the Bianchi identities of either G ' a V or 
Q^ v together with © or ©. 

Eqs. pTl [T2l [TB")) were the equations selected in Q for 
evaluating the numerical solutions, which is the simplest 
choice in the isotropic picture. Nevertheless, for plotting 
anisotropic cosmologies one does not have the possibility 
of selecting only the first order equations. To prepare to 
the next section, we will use the second order equations, 



2H + 3H 2 - 2HH 



N 



A 



NXY N 2 p r 



(14) 



a 2 y 6 



Besides p r and pb, two parameters, a and A, and two 
initial conditions, Yq and Hq, are still free. The proper 
way of fixing them is to search for the best fit directly 
from cosmological data fltl [l7l Il8j . Here we follow the 
original proposal and assume that this best fit should 
lead to values of p r , pb, A and Hq close to the ones found 
from the ACDM best fit. Any choice of these parameters 
close to the values found from WMAP5+BAO+SN [1] is 
perfectly consistent with this paper proposal, and we will 
take precisely those values — except for radiation, for it a 
mean value is adopted. Regarding the values of Yq and a, 
a is set to 0.99 (it should be close to one to allow solutions 
that are close to the standard cosmological model [|[), 
and Yq is set suc h th at the Hubble parameter directly 
evolves toward yfk/Z, leading to Y Q = 1.083 (H. Fig. 
(fTJ) shows a plot for the evolution of all the isotropic 
Eddington-Born- Infold functions. And Fig.Q shows the 
evolution of H, Hy and the pressure induced in (TT4"|) by 
the new fluid, it will be useful for the next section. 




2Hy + SHy — 2HyHx = 

A(l -a)X 2 f 2a 
1 — a 



1 

X 2 



a 

Y 2 



(15) 



together with (|13|) as the equations to be solved numeri- 
cally. With H x = X/X and H N = N/N. Eqs. ([TT]|T2)) 
constrains the initial conditions. 



FIG. 1: Plot for all the Eddington-Born-Infeld functions, with 
fl b = 0.05, Q r = 8.49 x 10" 5 , Q A = 0.72, H = 0.98, t = 1 
and q = 0.99. The values of Yb, Yo and Xo are inferred 
from the latter. The big bang happens at t = —0.787.... This 
deviation from zero is just outside the one sigma level error for 
the universe age Q. It is a negligible deviation, in particular 
since this model parameters have yet to be properly fit directly 
from the cosmological data 



nave yet to 



In order to plot the numerical solutions, we fix the 
temporal gauge with N = I and search for a consistent 
set of initial conditions. We will use Xg, Hq, Yq, Yq to 
denote the values of X, H, Y , Y at present time to, and 
we set a(to) = 1. The initial conditions will be set at 
t = to. Letting (jTTj) to set the value of Xo as a function 
of Hq and Yq, the Y initial condition reads 



Yo= ± 



At Y 4 



V6{- 



3H 2 



+ * + 



2a 



1 — a 



Y 2 



III. ANISOTROPIC PERTURBATIONS 

The approach of this section is to study the behavior 
of the anisotropies assuming that they are fluctuations 
around the isotropic solution. The isotropic quantities, 
whose analytical solutions are unknown, are analyzed 
numerically, while the anisotropic quantities are ana- 
lyzed analytically. In this model, as it will be shown in 
the following section, there are six initial conditions that 
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FIG. 2: Evolution of H = a/a, H Y = Y/Y and P/k = 
—AXY/a, considering the same parameters of Fig.{l]). All 
these three quantities are in above considered as adimensional 
quantities (using to — 1). Although the ratio of the pressure 
P by the energy density goes to zero as t approaches the big 
bang the pressure P does not go to zero. 



iliary tensor field 



(^)=diag(-A 2 Y? Y 2 Y 2 



(21) 



The equations of motion [5]) can be written in the 
form of (|18[ I19p with the appropriate metric variables, 
energy density and pressures. The energy density and 
pressures induced by the field on the equations for 
the g^i, dynamics read, henceforth fixing N — 1, 



Pq 



Pq 



fcA 



Y 1 Y 2 Y 3 



Xaia 2 a 3 



= -fcA 



XY 1 Y 2 Y 3 a* 
a\a 2 a 3 Y? 



(22) 



(23) 



Analogously, the metric induces a fluid on the q^ equa- 
tions whose energy density and pressures are 



fcA(l 



2a 




(24) 



parametrize the evolution of the anisotropies, therefore a 
pure numerical analysis for some sets of initial conditions 
does not unveil significantly the general picture. In 
this section it will be shown that, independently on 
the anisotropic initial conditions, the isotropic solution 
acts as an atractor for both early and late times, and 
that the Eddington-Born-Infcld isotropization has a 
characteristic pattern. 



The Bianchi I cosmology follows from general relativity 
with the homogeneous metric [l2| 



fcA(l - a) 



(9»v) = diag ( ~N 2 a\ a\ a\ ) 



(17) 



Introducing Hi = aJdi, its four equations of motion 
can be expressed as [231 



N 

H\H 2 + H 2 H 3 + HiH 3 — —— p, 

fc 



3 = 1 



HiH 2 H 3 
~ Hi 



(18) 



N 2 

-^■(19) 



These generalize the usual Friedmann equations for a flat 
anisotropic universe without rotation, with energy den- 
sity p and pressures pi . 

Using Ti. = ^ • Hj , the shears evolution is found from 

d N 2 

— {Hi - H k ) = (H N - m) [Hi -H k ) + — ( Pi - Pk ) . 

dt k 

(20) 



The cosmology proposed in @] is generalized to a 
Bianchi type I one by using the metric \Y1\ and the aux- 



2a 
1 - a 



1 

X 2 



3 



a 2 a 2 



/ j yi ^Y 2 

3=1 3 



(25) 



In the above we could have introduced another constant 
different from fc, but that is superfluous since the addi- 
tional C p a p connections do not couple to anything aside 
from gravity. The above energies and pressures reduce 
to their corresponding isotropic quantities if a, = a and 
Yi =Y (see El [II [13 [[5]). 

With the purpose of separating the isotropic contribu- 
tion from the anisotropic ones, we set 



Hi =H + e t 
H Yl =H Y +% 



(26) 
(27) 
(28) 



After the gauge is fixed (N = 1) a residual symme- 
try is still present, namely the rescalings of the spacial 
coordinates by constant factors, leading to the follow- 
ing transformations for the scale factors: di — > Cj 
and Yi — > Ci Yi, where the Cj's are arbitrary non-null 
constants. One should take care when dealing with 
isotropization directly from the metric components; in 
particular, that all the metric components approach the 
same value is not sufficient (nor necessary) for isotropiza- 
tion to happen; it can be just an artifact of the selected 
coordinate system. On the other hand, the quantities Hi , 
Hyi and p q are all invariant under the above symmetry, 
and the six independent quantities p qi — p q ., Hi — Hj, 
Hyi — Hy j fully characterize the anisotropies. 

To first order on the anisotropic perturbations 
(si, £j, 74), the equations of motion satisfied by Ji, Jiy and 
V are the same of their corresponding isotropic quanti- 
ties (H, Hy, p q ) if we set = Y<i s i = E,7» = 0. 
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This result only depends on the assumption that e, , 7, 
are small in regard to H, Jiy respectively. That is, &i is 
assumed to be small in regard to H, but its relation to 
other quantities is arbitrary. 

The evolution of the pressure anisotropies (anisotropic 
stresses) read 

j t (Pm ~ p«) = ^(6 - 6) = («y - W + H x ) te - 6) + 
+ 2V(e t - 7j - £j + lj ) + 2&(e< - 7 i) - 2&( £j - 7i ). (29) 

Whenever the anisotropies are small in regard to their 
corresponding isotropic values, the two last terms from 
(129")) can be dropped. 

All the dynamics of the anisotropies can be found from 
the two versions of ([20)1 (one for g^ v and the other for q^) 
and ([2"9")1 . In a standard Bianchi I cosmology there are two 
groups of two interacting anisotropic quantities, which 
are the shears (Hi — HA and the pressure anisotropics 
(Pi ~ Pj)- I n the Eddington-Born-Infcld extension to a 
Bianchi I scenario, the evolution of the latter depends on 
the physical shears (Hi — Hj) and the shears from q^ v 
(HYi — Hy j). There is not an independent second set of 
pressure anisotropies, see ([521 . 



A. Early time behavior 

To some extent, how the shears in the Eddington-Born- 
Infeld model will behave just after inflation depends on 
the role of the connections C a ^ during inflation. How- 
ever, at early times, if the anisotropic quantities are 
smaller than their isotropic counterparts, or the shears 
will decay or their derivatives will approach zero. There 
is no guarantee that the shears will decay from the early 
universe to the present time, and Figs. ([H [7]) illustrate 
that, but the isotropic solution acts as an attractor. 

For the shears modulus to increase it is necessary that 
1 6 — 61 > \SkTL(Ei — £j)\. The assumption of small 
anisotropies does not rule out this possibility at early 
times, for the pressure exerted by the unified dark fluid 
starts at minus infinity. Moreover, no upper limit in 
16 — 61 can avoid that inequality, since there is no lower 
bound for the shears. Nevertheless, if that inequality 
holds at early times, then 

4(6-6) ~ {n Y -n+H x m-^)+2V(- li + lk ). (30) 



(matter phase) or aT 2 (radiation phase) while the second 
to a -1 . 

Therefore, if after inflation the pressure anisotropies 
have a relevant contribution to the shears evolution, 
the former will quickly decay if the last term of ([30]) 
can be neglected. It might happen that the pressure 
anisotropies will not decrease just after inflation, for 
this to happen it is necessary that 2\V(— 7$ + 7 fc)| > 
I (Hy — H + Hx ) — £j ) | • Which implies, at early times, 
that \(H X - 3H Y )(% - 7fc)l » l A ( x ~ a )/p« (6 - 6)1, 
and hence, from ([201 13"2"|) , that (7, — jf-\ will decay. 

In conclusion, it is possible that the modulus of the 
physical shears increase at early times, but if this hap- 
pens the pressure anisotropies will decrease, reducing the 
shears derivative. The isotropic solution acts as an at- 
tractor at early times. 



B. Late time behavior 

For the late time behavior, another approach is capable 
of unveiling more details. Using the following relations 
that come directly from the definitions, 

X 2 Pq 

Pn - Pij = - fcA(1 _ a) (P 9i - P 9j ), (32) 

V = (H Y -H + H x )V, (33) 
p q = (3(H Y -H)-H x )p q , (34) 

it is possible to write the second order expressions for 
the evolution of the shears and pressure anisotropies as 
follows, 

d 2 , ( V\ , , d , 

-^(Si - Ej) = \-m + 2- J (Si - Ej) - 3H-( £l - £,) 

+ \(H Y -H + H x )(£i - fc) - 2^( 7l - 7,), (35) 
^(6 - 6) = (Hy - H + H X + 2 f - 2^^) (6 ~ 6) 
+ (H Y -H + H x )^i - Q + 2V(-H Y + W)(7< - 7i) 



The latter comes from 

\(H Y -H + H x m-W\> (31) 

\(H Y -H + H x )3kH(£ l - Ej)\ » 2\V(e t - Ej)\, 

since, as t goes to zero, H increases faster than \V\ (see 
also Fig. ^)) — the first evolves proportionally to a -3 / 2 



+2V(H Y - 4H + Hx)(Ei - Sj), 

- 7i) = [Hx ~ m Y + 2^—^ J (7i - 7,) 
+ (H X m Y A(n - 7i) - 2 VHl ' a) (E l - Ej ) 
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L(n Y -n-H x )(^-Q. (37) 

Pi 



Albeit less compact, the above form is useful for inter- 
preting the evolution of anisotropics as the interaction of 
three parametric damped oscillators. 

The differential equation x — —bx — lo 2 x describes a 
damped oscillator with damping factor b > and reso- 
nance frequency u>. The effective frequency of oscillation 
is i-\/4w 2 — b 2 . If this value is complex the x evolution 
is given by hyperbolic trigonometric functions, otherwise 
oscillations appear for x. We will use these observations 
as approximations to the late time behavior of (|35l36l37[k 
firstly without interaction, i.e. neglecting the two last 
terms in each of the above equations [30| . 

From the numerical solutions of the three damping co- 
efficients, b E — 3H, b 1 = —(Hx — 3Hy) and b^ — — (Hy — 
H+Hx), it is found that they are always positive, leading 
to a tendency towards isotropization. In Fig. the val- 
ues of b 2 ,b 2 and b 2 are plotted and compared to Au> 2 = 

4 (sH - 2f) , 4^ = 4 (-H x + m Y - 2^^) and 
4wf = 4 (-Hy+H- 

The effective frequency of oscillation i\/4a; 2 — b 2 for 
the physical shears and the pressures anisotropics be- 
comes real after t ~ 0.8 t and t ~ 0.2 i n respectively. 
These quantities oscillate with approximate asymptotic 
effective (angular) frequencies of 2/to and 3/(2io), and 
hence their half periods are 3/2 to and 2 to respectively. 
These late time low frequency oscillations can reduce 
dramatically the decrement of the physical shears, when 
compared to the standard Bianchi I exponential decay, 
but do not spoil the isotropization tendency, since the 
damping factors are always positive. 

The shears from q^ have an intrinsic oscillatory ten- 
dency at high redshift, but with low and quickly decreas- 
ing frequency. While for late times, it simply decay with 
an exponential-like behavior. 

When the two last terms of (|35[l3l)ll3"Tj) are considered, 
turning on the interaction among the three damped para- 
metric oscillators, the picture at t < 1 can change con- 
siderably depending on the initial conditions (see next 
section). During these times, due to the interactions, 
some oscillations with ill defined large period (> 0.5) 
can appear. High frequency oscillations cannot appear 
since all of the oscillators have low frequency. For latter 
times, the behavior is less dependent on the initial con- 
ditions. The physical shears, in particular, oscillate with 
a well defined low frequency and with decreasing ampli- 
tude, interacting only with the shears from q^ v (since 
Hy — H + Hx — > 0), whose natural tendency is to decay 
in an exponential way. 



20 



10 







-10 



V Z- 













! 1 
! 1 
1 / 


b e 2 l^ 


4cl> £ 2 /h 2 - 


1 1 
1 1 
1 1 


4oj e 2 /ff„ 2 


by 2 /iV - 


1 1 
1 1 
1 1 


bljm 




4w y 2 /s„ ! - 


1 .I. . « 


" 


0.0 0.5 


1.0 
t 


1.5 2.0 



FIG. 3: Evolution of the damping and resonance parame- 
ters on time. Generated from the same isotropic solution 
displayed in Fig.{T}. In particular, even though at late times 
this model isotropic version behaves as a standard de Sitter 
universe, its Bianchi I version does not display the typical ex- 
ponential decay for the shears, instead they decay oscillating 
(see Conclusions for other details). 



IV. NUMERICAL SOLUTIONS FOR THE 
EVOLUTION OF ANISOTROPIES 

In this section we present an approach to the numeri- 
cal analysis of this problem and numerical solutions that 
complement Sec. Ill results. For late times (t > 1) the 
anisotropics always decay oscillating with well defined 
low frequency, but many possibilities are left for the de- 
tails on the anisotropies evolution at both high and low 
redshift. This section displays some peculiar possibilities 
of this model, considering small anisotropies that fluctu- 
ate around the isotropic solution presented in Sec. II. 

The Eddington-Born-Infeld model with a Bianchi I 
background is described by seven equations of motion 
and seven unknowns, once the temporal gauge is fixed. 
Analogously to the isotropic case as presented in Sec. 
II, the equations selected to be solved numerically are 
the six second order ones (i.e., (fTT))) with the appropriate 
variables) and energy momentum conservation, 

We use X , HoijYoi and Yq ( to denote the values of 
X, Hi, Yi and at the present time t , and set dj(to) = 1- 
The first order equation (|I8[) for the metric determines 
Xq as a function of iJoi, Ya it while the first order equation 
for determines one of the Fch's a function of H^i-, Yoi 
and the two others Yq^s. 

In general, the value of Xq in the anisotropic case is 
different from the isotropic one. It remains to be imposed 
that the isotropic quantities Hq and Yq are the mean 
values of the anisotropic quantities. To this end, one sets 
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#o = IpLi-tfoi and F 3 = Y Q1 Y 02 Y Q3 .^ The latter, for 
small anisotropies, are sufficient for finding a X that is 
equal to the isotropic one and Hyo = | Yli=i H-Yoi- 

Figs. ([H [5]) exemplify the evolution of anisotropies at 
late times for two sets of initial conditions set at t$. Some 
details on the evolution at low rcdshift depend strongly 
on the initial conditions (e.g., the shears may or may not 
become null at t < 1), but at latter times their behavior 
follows the general pattern described in Sec. III. 
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FIG. 4: Evolution of all anisotropic quantities at late time. 
All the quantities in this plot refer to the corresponding 
adimensional ratios, i.e. {Hi — Hj)/7i, {Hy% — Hy j) /TCy , 
(Pn ~Vij)/'P- Only the two shears with greater magnitude 
are displayed, together with the pressure anisotropies and g M „ 
shears of same indices. All the anisotropies were set to be 
small at t = 1, but with a considerable magnitude difference. 
The physical shears were set to be 10 -5 , being close to the 
usual bounds, and the others were set to 10 _1 . This solution 
only has small anisotropies at late times. 
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FIG. 5: Like in Fig.Q, but all anisotropies are set to have 
about the same order ~ 10 -5 at to, and the sign of the pres- 
sure anisotropies at to was inverted, preventing the shears to 
continue to decay before reaching zero. 

From the numerical solution of the isotropic case (or 
from the anisotropic equations solution with isotropic 



initial conditions), one can infer the value of all model 
parameters at any time, and hence generate numerical 
solutions from initial conditions at any time. Figs. (0 
[7J show two sets of solutions for the shears evolution 
with nonstandard features. Both were generated from 
initial conditions set at matter-radiation decoupling time 
— which is a natural choice since the strongest shear ob- 
servational constraints come from the cosmic microwave 
background (CMB) [l2| ■ The dashed solution shows that, 
in spite of the isotropic solution being an attractor, it is 
possible that the shears have been increasing since the 
last scattering surface formation and up to the present 
time. This idea was explored recently in other models 
as a way for solving certain cosmological anomalies, e.g. 
[22L I24L HE]. The solid shears solution shows the pos- 
sibility that the shears could have been high at at the 
last scattering an high at present time, but with inverted 
signs. This possibility can significantly reduce the stan- 
dard observational constraints on the shears — at the 
price of adding one coincidence problem in cosmology, 
since the CMB will seem more isotropic at current time 
than at the past or the future. Note that all solutions 
of this model have shears that invert sign, since at late 
times the shears approach zero like damped oscillators. 
The special feature of this solution (the solid one) is when 
the first sign flip happens. 




0.0000 0.0002 0.0004 0.0006 0.0008 0.0010 

t 

FIG. 6: Evolution of two sets of shears with initial conditions 
set at decoupling time (about 10 to after the big bang). In 
above the instant t — was shifted to be the decoupling 
time and the shears are measured in Ho unities. For the 
dashed solution it was used, at decoupling: Hi — Hj ~ 10 -9 7i, 
H Yi - H Yj ~ 10~ 2 Hy, p qi -p qj ~ 10" 3 P. For the solid 
solution, the same as above, but Hi — Hj ~ 1Q~ 5 TL. 



V. CONCLUSIONS 

The Eddington-Born-Infeld (EBI) model Q is a new 
candidate for an unified dark sector and uses additional 
connections, instead of scalar fields, to model dark energy 
and dark matter. The approach with additional connec- 
tions makes its extension to Bianchi cosmologies nontriv- 
ial since it induces anisotropic pressures (stresses). 
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Flipping shears solulion 
Increasing shears solution 



FIG. 7: The late time behavior of the same solutions of 
Fig.©. 



In this work the Eddington-Born-Infeld cosmology was 
extended to a Bianchi type I cosmology and it was shown 
that i) the full dynamics of anisotropics depends on 
two noninteracting groups of three anisotropic quantities 
(the physical shears, the shears from q^ U: and the pres- 
sure anisotropies); ii) for both early and late times the 
isotropic solution acts as an attractor, implying that the 
shears tend to decay, and that observational constraints 
on the shears can be fulfilled [HI, Hi) the shears 
do not decay monotonically, but they typically oscillate 
with decreasing amplitude; iv) letting some fine tun- 
ing on the initial conditions that parametrize this model 
anisotropies, we could construct two solutions qualita- 
tively different from usual anisotropic models Figs. (0 
[7]); specifically, these are capable of avoiding standard 
theoretical bounds on the shears. 

Since here we have evaluated anisotropic fluctuations 
around the isotropic EBI solution found in Q, some of 
our results rely on the same assumption considered in 
[U, namely that some of the EBI isotropic parameters 
should have values similar to the corresponding ones from 
ACDM 4]. We stress that further EBI developments 
may lead to different parameters, for the Ref. [4j results 
were found assuming ACDM. On the other hand, many 
qualitative aspects of the Bianchi I EBI cosmology here 
presented do not depend on the precise values of those 
isotropic parameters; in particular, from Eqs. (|2QI29[) or 
(|35I36I37[) one sees the shears oscillations is a general fea- 
ture of this cosmology 3l| . There is yet much to explore 
in anisotropic EBI systems, we should return to them in 
a future work. 

Collins and Hawking [2(| analyzed the issue of stability 
of the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
cosmology under anisotropic perturbations, and found 



that the universe is unstable at late times in some Bianchi 
models. However, a cosmological constant at late times 
was not considered in their anal ysis , which is known to 
be sufficient to lead to stability [21(. At late times, the 
isotropic EBI model effectively behaves like FLRW with 
cosmological constant. Nevertheless, the anisotropic per- 
turbation induces a modification in the effective cosmo- 
logical constant and it starts to exert anisotropic pres- 
sure. It can be inferred from Eq. ([29)) and Fig. ([3]) that, 
if a shear perturbation occurs at late time, the shears 
will induce anisotropic pressure which will prevent the 
exponential decay [32j, but the shears will decay follow- 
ing a damped oscillatory behavior; hence EBI cosmology 
is also stable at late times. 

Recently anisotropic forms of dark energy have been 
receiving increasing attention 

d3[H0[2ij& These 

proposals, besides constraining what dark energy can and 
cannot be, are generating their first observational conse- 
quences and perhaps are the explanation behind some 
anomalies found in the standard cosmological model. 
The Eddington-Born-Infcld model naturally generates an 
anisotropic form of dark energy that is considerably in- 
tricate. It is not the purpose of this paper to pinpoint 
a set of initial conditions capable of solving a particu- 
lar anomaly; but we would like to stress that standard 
bounds on the shears need not to be satisfied in this 
model. 

The present work is a step toward the proof that the 
EBI model is in conformity with observational data for 
modeling dark energy and dark matter. Of special rele- 
vance to this proof is the analysis of density fluctuations 
and consequences for the CMB anisotropies 16], since 
this step the Chaplygin gas approach could not overcome 
as a unified dark sector model [7[. If the EBI model can 
indeed be a realistic "quartessence" model it has yet to 
be proven, but this paper results do not depend on that. 

Considering the original motivation Q , the EBI model 
should be more than a phenomenological model for dark 
energy and dark matter; the role of the C^ p connections 
during inflation have yet to be evaluated and other the- 
oretical issues are being evaluated. 
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